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ABSTRACT 

We present three transits of GJ 1214b, observed as part of the Apache Point 
Observatory Survey of Transit Lightcurves of Exoplanets (APOSTLE). We used 
APOSTLE r-band hghtcurves in conjunction with previously gathered data of 
GJ 1214b to re-derive system parameters. By using parameters such as transit 
duration and ingress/egress length we are able to reduce the degeneracies be- 
tween parameters in the fitted transit model, which is a preferred condition for 
Markov Chain Monte Carlo techniques typically used to quantify uncertainties in 
measured parameters. The joint analysis of this multi-wavelength dataset con- 
firms earlier estimates of system parameters including planetary orbital period, 
the planet-to-star radius ratio and stellar density. Estimating the absolute mass 
and radius of the planet directly depend on how various stellar parameters are 
derived. We fit the photometric spectral-energy distribution of GJ 1214 to derive 
stellar luminosity and estimated its absolute mass and radius from known mass- 
luminosity relations for low-mass stars. From these derived stellar properties and 
previously published radial velocity data we were able to refine estimates of the 
absolute parameters for the planet GJ 1214b. Transit times derived from our 
study show no evidence for strong transit timing variations. Some lightcurves we 
present show features that we believe are due to stellar activity. During the first 
night we observed a 0.8% rise in the out-of-eclipse flux of the host star lasting 
approximately 3 minutes. The trend has a characteristic fast-rise exponential 
decay shape commonly associated with stellar flares. During the second night we 
observed a minor brightening during the transit. Due to its symmetric shape we 
believe this feature might have been caused by the planet obscuring a star-spot 
on the stellar disk. 
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Introduction 



Exoplanets detected using both transits and the radial velocity (RV) technique offer the 
unique opportunity to measure many of their physical properties. Transits can place limits on 
the radius of an exoplanet given the radius of the star, and transits can also constrain orbital 
parameters such as inclination and orbital period. When coupled with RV measurements, 
we can constrain the mass (Mp), and the average density of a planet (pp). We may therefore 
probe t he interiors of transi t ing exoplanets and constrain bulk corn position and formation 
models (IFortney et al.l l2008l : [Torres et al.l l2008t iBaraffe et al.l 120081). Here we present a nd 
interpret three lightcurves of the transiting planet GJ 1214b (jCharbonneau et al.ll2009l ). a 
planet unlike any in our Solar System. 

Prior to 2009, all transiting exoplanets were found with radii consistent with the giant 
planets in our Solar System, i.e. they had gaseous envelopes. As surveys improved, they 
began to reach sensitivities which could detect smaller rocky planets. In our Solar System, 
Neptune has a mass of 17 while the Earth is the largest terrestrial planet, suggesting 
that the transition mass between rocky and gaseous planets lies somewhere between these 
values. With no such "transition" object in the Solar System, w e must rely on theory, which 
predicts that ~ 10 is th e critical mass (IPoUack et al.lll996l). althoug h it could be as low 



as 2 Mq^ fllkoma et al.l 120011 ) or as high as 16 Mg. (]Lissauer et al.l 120091 1. The 10 M® limit 
should therefore be seen as a sort of median of theoretical results, and not a true boundary 
between rocky and gaseous worlds. The discovery of transiting planets between 1 and 17 
M^therefore provides critical insight into the planet formation process. 



Two transiting plane ts are now know n with lMm< Mp < 10 M^ ; CoRoT-7b fjLeger et al. 



20091 : iQue 



oz et al. 



20091 



]), a nd GJ 1214b Charbonneau et al. (2009). The former appears to 



be rocky ( Leger et al.l 2009), while the la tter may contain significant amounts of water or a 



gaseous envelope (IRogers &: Seageiil2010l ). In order to address the ambiguities in the planet 
formation process, the radii of these planets must be measured accurately and precisely. The 
CoRoT satellite already surpasses all other proj ects in the ability to make follow-up obser- 
vations of CoRoT-7b (153 transits reported by iLeger et al.l (120091 )). The relatively recent 
detection and hence the dearth of follow-up measurements on GJ 1214b l ed us to focus on this 



objec t; a Super-Earth planet orbiting M dwarf star 13 pc from E arth (IChar 



20091 ) . The planetary nature of this transit has been confirmed by ISada et al 



jonneau et al. 



mm 



The discovery of such a small-sized planet bodes well for transit searches for habitable 
planets around M-Dwa rfs. The low luminos i ty of M-Dwarfs rn ean their habitable zones are 
very close t o the star (iKasting et al.l Il993l : ISelsis et al.l 120071 ). This increases the transit 
probability (iBorucki fc Summers! 119841 ). and hence some of the first planets to be character- 
ized as rocky and in the habitable zone may well be transiting planets around M-Dwarfs. 
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Although GJ 1214b most hkely possesses a Hydrogen-rich envelope, the detection of small 
planets such as this heralds the discovery of rocky habitable worlds if they e xist. M-Dwarfs 
make up a very large fraction of the stellar component of the Milky Way (IMiller &: Scalo 



19791 : iReid et al.ll2002l ). so the prospect of habitable planets around M-Dwarfs raises the very 



interesting possibility that life-bearing planets may be fairly common in the galaxy. 

In this paper we report observations of three transits of G J 1214b in 2010. In § [2] we 
outline our observations and data reduction techniques; in § E] we describe our lightcurve 
model and the use of Markov Chain Monte Carlo (MCMC) techniques to constrain system 
parameters from single and multi- wavelength data. In § |4] we report the absence of strong 
transit timing variations (TTV). In § [5] we describe the derivation of various stellar and 
planetary parameters for the GJ 1214 system. In § Elwe discuss the detection of a flare and 
a possible spot-crossing event. Finally, in §[7] we summarize our findings. 



2. Data 

2.1. APOSTLE Observations 

We observed three transits of GJ 1214b on UT dates 2010-04-21, 2010-06-06 and 2010- 
07-06 using the ARC 3.5m Telescope at Apache Point, New Mexico. All observations were 
made with Agile, a high-sp eed time-series CCD photometer based on the design of Argos 



( iNather &: Mukadanj|2004j ). Agile is a charge transfer CCD, that collects photons from the 
target at a 100% duty cycle. During transit observations the charge on Agile was read out 
at 45 sec intervals using GPS-synchronized pulses with an absolute timing accuracy of less 
than a millisecond. Observations were made using Agile^s medium gain, slow readout mode 
with frames binned by a factor of 2, yielding a plate-scale of 0.258 "/pixel. All APOSTLE 
lightcurves pr esented in this work were collected using the r-filter, which is similar to the 



SDSS r filter (iFukudta et al.lll996h . 



During observations we defocused the telescope to spread the stellar point-spread func- 
tion (PSF) across multiple pixels in order to minimize systematics introduced by improper 
flat-fielding and to increase the integrated counts from the star. Count rate stability is af- 
fected by variable conditions during observing. We monitored the count rate and adjusted 
the telescope focus by small increments to raise or lower the maximum counts on the bright- 
est star. By adjusting the focus we kept the maximum counts between 40k and 55k ADU. 
We did this because the instrument is known to show a non-linear response at count-levels 
greater than ~55k, and saturates at 61k counts. During stable conditions these adjustments 
were made every 10 to 20 minutes while during poor conditions the adjustments were more 
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frequent, at roughly 5 minute intervals. Less than 7% of frames on a given night were lost 
to saturation. For all our observations the comparison star USNO-Bl.O 0949- 0280051 was 



the b rightest object within Agile^s field-of-view, 95" to the north of GJ 1214 (IMonet et al. 



20031 ). This object was brighter than GJ 1214 by a factor of 1.4 in the r-band. It was also 



the only star brighter than GJ 1214 in Agile^s field-of-view, and hence was the only one 
used for differential photometry. On each night we also collected twilight sky flats and dark 
frames (which were also used to correct for bias). 



2.2. Image Reductions 

We used a customized data reduction pipeline, written in Interactive Data Language 
(IDL) to process Agile data. It performs standard image processing steps like dark subtrac- 
tion and fiat fielding, but also implements non-linearity corrections unique to Agile. The 
pipeline also creates an uncertainty map of the processed images by propagating pixel-to- 
pixel errors through each step of the reduction. In addition to the Poisson photon counting 
errors and read-noise from the science images, the pipeline propagates the variance on the 
master dark and master flat during the reduction. Errors were also propagated for those 
pixels where the counts exceeded the non-linearity threshold (~55k counts) using the un- 
certainties in the empirically derived non-linearity correction function. Typically 1-5% of 
frames went above this threshold. The correction factor for counts on non-linear pixels was 
typically 1.45% or less. 



2.3. Photometry 



We used SExtractor (IBertin fc Arnoutslll996l ) to derive initial centroids of our defocused 



stars. Coordinates obtained from SExtractor were then used for circular aperture photometry 
with the PHOT task in IRAF's NOAO.DIGIPHOT.APPHOT package. We derived flux 
estimates from a range of circular apertures with radii between 5-50 pixels, at intervals 
of 1 pixel. An outlier-rejected global median was used as the sky estimate. An optimal 
aperture was selected where the RMS in the out-of-eclipse lightcurve was minimized. For 
GJ 1214 data this aperture was typically 15-16 pixels in radius. This size was roughly 
4 times the Half- Width at Half Maximum (HWHM) of the PSFs. To derive photometric 
errors we extracted counts from the error frames using the same centroids and apertures 
used for photometry on the target frames. Estimating photometric errors in this manner 
yielded uncertainties that were greater (by 50%) than the default errors reported by PHOT. 
It is useful to note that this method takes into account sources of error which are otherwise 
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ignored by standard photometric techniques, e.g. fluctuations in the flat field and dark 
frame, and is thus more thorough. Photometric precision for our GJ1214 lightcurves was 
typically ~ 0.001 magnitudes from 45sec exposures. 



2.4. Time Coordinates 



The lAU recommended time co ordinate for t he proper analysis of event timing data is 
TCB (Barycentric Coordinate Time. lKaplarul2005l ) . There are two primary advantages of us- 
ing this time coordinate: 1) TCB is the best approximation we have to taking measurements 
from an inertial reference frame. In this system, the time-of-arrival of photons are described 
as measured by a clock at the solar system barycenter (SSB) assumin g gravitational effect s 
on light travel time due the largest solar system bodies are rem oved (lEdwards et al.l 120061 ). 



2) TCB is based on the SI definition of the time unit - second ( IStandish 



19981 



W e followed th e presc riptio n described in ISeide 



mann fc Fukushimal ( 1l992l ). IStandish 



( 119981 ) ■ iHobbs et al.l ( 120061 ) . and Edwards et al.l ( 120061 ) to transform the mid-exposure UTC 
time (Coordinated Universal Time) into TCB: 



TCB{JD) = UTC{JD) + Ccor + Gcor + 



:ii 



where TCB (JD) and UTC ( JD) are the Julian Day representations of TCB and UTC respec- 
tively; CcoT is the clock correction term (used to convert UTC to Terrestrial Time (TT)), Gcor 
is the geometric light travel time cor rection (or Roemer delay), and Ernr is relativistic cor- 
rectio n (the Einstein delay) (see also, iHobbs et al.1l2006[ lEdwards et al.ll2006l : lEastman et al.l 
2010). The method is also desc ribed in detail by Eastman et al. ( 201ol ) with reference to 
exoplanetary transits. However, lEastman et al.l (120101 ) recommend the closely related TDB 
(Barycentric Dynamical Time) system instead of TCB for reporting ti ming measurern ents. 
The difference between these two coordinates is very slight. However, IStandishI (119981 ) dis- 
cusses how the TDB system does not represent a physical time measure but is related to 
TCB by a constant offset and scaling factor. Due to this the time s represented by T DB are 
in units that are subtly different from the widely used SI second (IHobbs et al.l 120061 ). UTC 
and Terrestrial Times (TT) reported by observatory and GPS (Global Positioning System) 
clocks are typically in the SI system. 



2.5. Other lightcurves 



The MEarth team graciously shared the lightcurves presented in ICharbonneau et al. 

(120091 ) ■ which were either from the array of MEarth (8 x 40cm) telescopes or 2;-filter obser- 
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vations from the Fred Lawrence Whipple Observatory's 1.2m telescope (FLW01.2m). We 
included these data in our multi-wavelength analysis of system parameters. We modified 
the time coordinate to the TCBj/) system as described in §2.4^ and when data on nuisance 
parameters were available, we included them in our detrending analysis. 



3. System Parameters 
3.1. Transit Lightcurve Model and Detrending 



APOSTLE tran sit lightcurves of GJ 12 14b are shown in FigurelH We use the lightcurve 
models described in iMandel &: Agoll (120021 ) assuming quadratic stellar limb-darkening and 
a circular planetary orbit. Our transit code can fit for multiple transits simultaneously 
and allows for multiple sets of limb- darkening coefficients when data gathered using different 
filters are used. For the analysis presented in this paper we used models which fit parameters 
in either one of the following two sets: for lightcurves from a single-filter, we used Set 1: 
= {tTjtc, D,vi,V2,T0i^^^^^} and for the joint analysis of lightcurves gathered at different 
wavelengths, we used Set 2: 62 = {tT^tc, Rl/Rl,VlJ,,,Np,V2J...NF^T0i,,,NT}■ 
^he parameters tr and to are approxim ately the transit du ration and ingress/egress 
duration respectively (same as T and r from ICarter et al.ll2008l ). The transit duration is 
defined as the time between the middle of ingress to the middle of egress. The ingress 
duration describes the time between the start of the echpse till the planet has completely 
crossed the limb of the star. The duration of egress is assumed to be the same. In the limit 
of zero eccentricity orbits these two parameters can be approximated as. 



ti 



t 



V 

Rn 



G 



(2) 



(3) 



where the orbital speed v = 27ra/P and the impact parameter b = asini/R^,. The terms a, 
P, i, Rp and R^ are the s emi-major axis, or bital period, inclination, planetary radius and 
stellar radius respectively ( ICarter et al.ll2008l ). 



We assume quadratic limb-darkening as described in lMandel fc Agoll (120021 ). and assume 
the transit depth approximately changes as. 



Rl 



m = '-§{1 - nil - VI 



62) - 72(1 - va^)2) 



(4) 
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here 71 and 72 are the quadratic hmb-darkening coefficients from lMandel &: Agoll ( 120021 ) and 
Rp/Rl is the square of the planet-to-star radius ratio. Solving for the impact parameter in 
equations [2] and [3] gives b = a/1 — {It /tG){Rp/ Ri?)- So given It, to and the limb-darkening 
coefficients, the maximum depth at mid-transit (D) can be determined by solving a sextic 
equation in R^/R^. The difference between the two fit parameter sets (^1 and ^2) are the 
variables D and Rp/R^. Since the limb-darkening coefficients are dependent on the waveband 
used for observations, Rp/ Rl is commonly used to fit for the depth of transit lightcurves, but 
D might be better suited for constraining the transit-depth of single-ffiter data, especially 
when the star is expected to be strongly limb-darkened. 

The terms vi and f 2 are linear combinations of the quadratic limb- darkening coefficients; 
■^^1 = 71 + 72 and V2 = 71 — 72- These linear combinations were used since it is known 
that directly fitting for limb- darkeni ng coefficients resu lts in strongly anti-correlated error 



distributions for transit parameters ( iBrown et al.ll200ll ). In order to avoid unphysical limb- 



darkening proffies we applied the bounds vi + V2 > and < fi < 1. The TO terms are 
the times of transit center. The subscripts i...NT and j...Np are used to denote multiple 
transits (A^t) and multiple filters (Np) respectively. Together these parameters define the 
model transit lightcurve. 

There are many systematic trends which might be introduced over the course of ob- 
serving that cannot be accounted for using the reduction protocol described in § 12.21 For 
example, differential extinction due to airmass variation or photometric variation due to cen- 
troids wandering over pixels of varying sensitivities on an imperfectly flatfielded image. So 
for each image we extracted a set of nuisance parameters which were then used to compute 
a correction function {i.e. detrending function) to remove systematic trends. We modeled 
this function as a linear combination of nuisance parameters: 

Fcor,i = 0,0 + OfcXfc^j, (5) 

k=l 

where Xk^i are the nuisance parameters and ak^i are the corresponding coefficients. A typical 
set of nuisance parameters included (i) the airmass, (ii) the centroid positions of the target 
and reference stars, (iii) the local sky around the target and comparison stars, (iv) the 
global sky and (v) the total counts in the area that defined the photometric aperture on the 
masterdark and (vi) the master skyflat. 

For each night, the entire lightcurve is normalized to 1 by an initial best-ffi transit 
model. Then we fit for the coefficients of the correction functions using these residuals 
and a generalized linear least squares minimizer. The oq term is set to a constant, which 
in this case is 1 because of the way the lightcurve is normalized. We found that the most 
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significant systematic trends in the residuals were correlated with airmass and the counts 
on the masterflat. This tells us that differential extinction and imperfect flatfielding are the 
two greatest sources of systematic effects. 

Once the correction function and the model are derived we compute the goodness of fit 
to our data as: 

""^^ fO.,-MAe)-F,ori? 

I ' 

where the Oj and are the observed data and associated errors, Mj is the transit model and 
Fcor,i is the detrending function. For all subsequent optimization with either Markov Chains 
or a Non-Linear Minimizer, Equation |6] is used to evaluate the goodness of fit. One must 
note that we devised routines such that the correction function is recomputed along with 
the model lightcurve for each step in the Markov chain or each iteration in the Non-Linear 
Minimizer. 

In §[6] we note the possible evidence for stellar activity in some of our lightcurves. These 
features were seen in lightcurves after reductions and detrending, and we believe they are 
not associated with any systematic effects in the data. We excluded these points from our 
analysis of system parameters. 



3.2. Markov Chain Monte Carlo 



Bayesian inference techniques like Markov Chain Monte Carlo have become a popular 
tool for constraining system parameters from observational data. We used the Metropolis- 
Hastings (M-H) algorithm, a well known M CMC method, t o con strain the u ncerta inties 
for fitted parameters in our transit model. iTegmark et al.l (120041 ) and iFordI (120051 ) have 
very good descriptions of the algorithm and it s app li cation to relevant astronomical data. 
Our prescription is closest to that described by iFordI (120051 ). We approximate the posterior 
distribution and joint probability distribution of our model parameters given the observed 
data O, as P(6> | O) oc P{0)P{O \ 6) oc e-^^^^''/^. Our MCMC routines use the standard 
stepping and selection rules of the M-H algorithm. The jump functions for our parameters 
are of the form, 

e,+^ = e, + G{Q,<Tl)f (7) 

where and ere are the vectors of model parameters and their associated step-sizes respec- 
tively and G(0, cxg) is a random number drawn from a normal distribution with a mean of 
and a variance of a^. The factor / is an adaptive step-size controller which is used to guide 
the chain to the optimal acceptance rate. For the case where a jump is performed for the 
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entire v ector of model para meters, it has been shown that the optimal acceptance rate is 
~ 23% ( IGelman et al.ll2003l ). This desired rate is achieved by adjusting the step-size con- 
troller (/) every 100 accepted steps according to fmw = 434 fni.d/Ntrin i.'^, where Ntrini-: are th e 



numb er of steps attempted for the last 100 accepted steps (see also ICoUier Cameron et al. 



By varying the entire vector of model parameters and applying a single step-size modifier 
we risk the situation of using mismatched step-sizes and undersampling targeted posterior 
distributions. Statisticians have shown that well co nstructed chains wi ll properly sample 
posterior distributions at the correct acceptance rate (IGelman et al.ll2003[ ). However, the ac- 
ceptance rate is guided by the location of the chain in parameter space. The M-H algorithm's 
basis for sampling the posterior distribution lies in the fact that steps in low probability (large 
X^) regions of parameter space are accepted less often than those in high probability regions. 
For example, if we make a poor choice and select a starting step-size for one parameter to 
be too large relative to the other parameters in the ensemble, this parameter will traverse 
between regions of high and low probability faster than the rest. The acceptance rate in this 
case is biased by jumps in one parameter, and the chain will have poorly sampled posterior 
distributions for the remaining parameters. The key to a well-constructed chain is to choose 
the relative starting stepsizes for parameters such that they all roam high and low proba- 
bility regions of parameter space at roughly the same rate. To find this ideal set of starting 
step-sizes we ran a set of exploratory chains (40,000 iterations), stepping only one parameter 
at a time, until the step-size controller settled the simulati on to an acceptance rate of ~ 44% 
(the optimal acceptance rate for the one-dimensional case, IGelman et al.ll2003l ). Most of the 
chains reached this acceptance rate at an iteration between ~1000-7000. The values of the 
adaptive stepsize controller (/) near the end of these short runs gave us an idea of how much 
we under/over estimated the starting stepsize for a given parameter. Typically, these "final" 
/ values multiplied by the starting stepsize (for the short runs) proved to be good choices 
for the starting stepsizes to be used for the larger, multi-parameter MCMC runs. 

We ran 10 long MCMC chains for different combinations of parameter sets and lightcurves. 
Table [1] lists the names of the chains, the corresponding lightcurve data, the parameter sets 
used and some statistics from our post-run analysis of the chains. The chains are num- 
bered from 001 to 005 which represent the 5 different data sets used - listed in column 2. 
Chains 001, 004 and 005 are single-filter data sets corresponding to APOSTLE, MEarth and 
FLW01.2m observations respectively. These were used to test the models where the transit 
depth p was fit, for the h i ghly limb-darkened case (see § 13. ip . Chain 002 was simply a 
redo of ICharbonneau et al.l ( l2009l )'s analysis with our transit model and MCMC framework. 
In chain 003 we fit for transit parameters using all 3 data sets. The tags 'a' and 'b' de- 
note whether the limb-darkening coefficients were left fixed or open respectively. For chains 
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where the hmb-darkening parameters were fixed ('a') we either chose values from the hter- 
ature or used values which we re foun d to be suitable by others. For the APOSTLE r-band 



dataset we chose values from IClaretl (120041 ) for a 3000K star: {vi^apostle^V2^apostle) 



(0.908,0.305) ■ For the ME arth and FLW01.2m data we used unpublished values used by 



Charbonneau et al.l ( l2009l ) in their fit (priv. comm. P. Nutzman). For the MEarth "filter" 
we used {vi^m Earth, V2, MEarth) = (0.145,0.639) and for the FLW01.2 z-band filter we used 
{vi^FLWoi.2m,'V2,FLWoi.2m) = (0.404,-0.289). Since it is well known that these parameters are 
highly degenerate, the 'b' chains were run as a test of how well these lightcurves could be 
used to constrain stellar limb- darkening. The parameter set corresponding to each chain is 
listed in column 3 of Tabled For the chains with multi-wavelength data, like 002 and 003, 
the parameter set 02 was used, and for single filter data we used the set 9i (see § 13. ip . 

The typical number of iterations used for long MCMC chains ranged between 1-2.5x10®. 
These computations took a total of 80 CPU hours to complete on Linux workstations. For 
each chain, we selected only those steps where the acceptance rate remained roughly within 
5% of the optimal acceptance rate. Cropping the chains in this manner meant, many of the 
initial steps were discarded as most chains started with acceptance rates which were far from 
the optimal rate. 

The M-H algorithm's jump function is only dependent on the previous location of a step 
in the chain, so it is not unusual for sequential points in the chain to be correlated. The 
dimensionless autocorrelation function provides a good assessment of how many independent 
points there ar e in a chain. We computed this for each open parameter per chain, using the 



prescription of iTegmark et al.l (120041 ). and derived the correlation and effective lengths for 
each chain. These are presented in Table [H The effective length (the chain length divided 
by the correlation length) is a measure of the number of independent points in a chain and 
must be large (^ 1) for the errors derived from a chain to be meaningful. A large number 
of independent points signifies a statistically significant, well-sampled posterior distribution. 
We found that chains where the limb- darkening was fixed ('a' chains) had the greatest 
effective lengths and chains where limb-darkening parameters were left open ('b' chains) had 
the shortest effective lengths (see Tabled]). Chain 003b had the lowest effective length. This 
chain also had the largest number of open parameters (including 3 pairs of limb-darkening 
parameters). We believe the large number of limb-darkening parameters resulted in slow 
convergence for this chain. We discuss the significance of the resulting parameters and 
errors in the following section. 
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3.3. Parameters and Errors 



System parameters derived from our analysis are presented in Table |2] and Table |3] for 
the single and multi-wavelength datasets respectively. Directly fit model parameters are 
listed in the table sub-section titled 'Model'. These correspond to the variables described in 
section § 13.11 as part of the parameter sets 0^ and O o- The values liste d in the tables were 
obtained using the minimization package MINUIT (j James et al.lll994j ). The from the 
best-fit model and the degrees-of-freedom (DOF) are listed in the last two columns of Table 
[U The ensemble of points from the posterior distributions of 'Mode l' points were then used to 
compute posterior distri butions of various 'Derived' parameters (jSeager fc Mallen-Ornelas 



2003 



Carter et al.l 120081 ). We chose to present the following 7 derived quantities: i. the 
planet-to-star radius ratio ii. the orbital period, iii. the impact parameter (6), 

iv. semi-major axis in stellar radius units (a/i?^), v. orbital inclination {i), vi. orbital 
velocity (f) normalized by stellar radius, and vii. the stellar density (p*). The errors for 
both the 'Model' and 'Derived' parameters were then computed by sorting these data and 
choosing the 68.3% confidence intervals to represent the la uncertainties. For cases where 
the parameter values were distributed asymmetrically around the best-fit value (at a level 
> 10%), we present upper and lower uncertainty estimates. 

Figure |2] shows the joint-probability distributions of directly fit model parameters from 
chainOOla. The cross-hairs on each sub-plot mark the best-fit values. We can see that the 
distributions of 'Model' parameters shown in Figure |2] show little or no correlation between 
each other. The use of parameters like transit duration (tT-) and ingress/egress duration 
{to) in our lightcurve model has resulted in a parameter set with few degeneracies. Fig- 
ure E] shows results from the same chain, but of derived parameters. Here we can see that 
most of the derived parameters show strong correlations with each other. When character- 
izing lightcurves, it is fairly common practice to use some of the pararn e ters listed here as 
'Derived' as the directly fit model parameters (e.g. iHolman et al.l l2006l : IWinn et al.l 120071 : 



Collier Cameron et al.ll2007l ). Also common is the use of MCMC to obtain error estimates 



from such models. The presence of such correlations makes the interpretation of MCMC re- 
sults challenging. These degeneracies can result in 1) chains that have short effective lengths 
(slow convergence) or 2) incomplete sampling of posterior probability distributions. Both 
problems can be reduced by running longer chains. However this solution is not always 
practical as the number of required steps may be prohibitively large or unknown. In § 13.21 
we noted the problem of slow convergence for chains when the limb-darkening was allowed 
to vary. Similarly, the problem of insufficient sampling of posterior distributions is common 
when there are "banana-shaped" degeneracies between parameters. In Figure [3] we see such 
a degeneracy for the case of b vs. p^, b vs. a/R^, and b vs. Rp/R^,. In such situations a chain 
may get stuck in "valleys" of low and fail to sample other regions of parameter space. We 
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have shown how both these issues can be avoided by choosing a parameter set that is free of 
mutual degeneracies. Figure [2] graphically confirms this and column-6 (Corr Length) in Ta- 
bled] quantitatively establishes that chains with parameters that show no mutual correlations 
(chains 'a') converge quickly. 

Figure m shows our results from chainOOlb. This chain is similar to chainOOla, except 
the limb-darkening parameters Vi and f 2 were allowed to vary. As mentioned in § 13.11 we 
applied bounds {vi -1-^2 > and < f 1 < 1) to the limb-darkening parameters so as to avoid 
unphysical limb-darkening profiles. We note that correlations exist between the two limb- 
darkening parameters (FigH] panel - viapostle vs. V2apostle)- In addition, parameters 
that showed uncorrelated distributions in Figure [2l now show signs of being affected by 
degeneracies. The most strongly affected are the parameters t^, to and D. Not only do 
these parameters show degeneracies with the limb-darkening parameters, but they also show 
correlations between each other. The degeneracy between the ingress/egress duration {tc) 
and the limb- darkening parameters can be understood simply by the fact that both affect the 
overall shape of the lightcurve. The limb-darkening can change the profile of the ingress and 
egress regions of the lightcurve, while to determines the start and end points of ingress and 
egress. So given the error bars and scatter in our lightcurves, a range of limb-darkening and 
to values can produce good fits. Our inability to constrain these parameters simultaneously 
seems to suggest that milli-magnitude photometry might not be sufficient for the complete 
characterization of transit lightcurves with limb- darkening. These degeneracies also strongly 
affect the best-fit parameters derived from these chains. For example, in the list of results 
for 'b' chains in Tables |2] and |31 we find that values for tc and produce very discrepant 
estimates for derived system parameters (like b and p^). We conclude that results from the 
'b' chains are in general unreliable, so we elect to discuss results in which the limb-darkening 
parameters are held fixed. 



3.3.1. Which parameter set to consider? 

We analyzed a combination of five different datasets and two parameter options (chains 
'a' and 'b', see Tabled]). Out of the five 'a' chains three chains, 001a, 004a and 005a were 
single filter datasets (see Table [2]). These serve as simple checks on our transit model and 
MCMC framework. They also highlight the use of the maximum transit depth (D) instead 
of the square of planet-to-star radius ratio {Rp/ Rl) for single filter data. We find that D has 
weaker covariance with to and than does R^/Rl- Due to the small number of lightcurves 
for each of these three sets we see, not surprisingly, that parameters such as the planet-to- 
star radius ratio {Rp/R^) and period are constrained to low precision. The chains numbered 
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002 were simply a redo of ICharbonneau et al.l (120091 ) 's analysis with ou r transit model and 



MCM C framework. We find the results are in good agreement with ICharbonneau et al. 



( 120091 ) 's findings. The largest parameter set presented in this work is chain003a. This set 
represents the joint analysis of multi-wavelength data u sing our transit model fram ework, 
and is in good agreement with the results presented in ICharbonneau et al.l (120091 ). From 
Table [1] we can safely say that the MCMC analysis for this chain is quite reliable. Thus the 
results from chainOOSa are the best to consider from this work and are used in the discussions 
that follow. 



4. Transit Timing Variations 



At the time of writing this paper, 10 tra nsits have been reported for GJ 1214b (^including 



Charbonneau et al 



2009; 



Sada et al.l l2010l ). Transit times reported in ICharbonneau et al. 



(l2009h and ISada et al.l (120101 ) were in HJOt/Tc*, i.e. JD representation of UTC corrected 



for the light travel time delay to the solar system barycenter. These were converted to ou r 
preferred BJDrc^ time coordinate. For the transits observed by ICharbonneau et al.l (120091 ). 
in Table [3] we present both the conversion of values cited by them and the transit times 
derived from our lightcurve fit to their data (see columns ' ICharbonneau et al.ll2009l ' and 
'chain003a'). We chose to use the times which resulted from our analysis (chain003a) since we 
fit for many transits simultaneously. These BJOt^b times are not far off from the converted 
transit times. The largest discrepancy we found was ~ 15 sec for the transit T0.T4. This 



differe nce is not surprising since it is the very first transit observed by ICharbonneau et al. 
(120091 ) and so, has the least precise measurement (a ~40 sec). With all available transit 
times we refit the ephemeris for GJ 1214b: 

TT{NTi) =Tzp + Px Nn (8) 

where TT is the expected transit time, Tzp is the zero-point of the transit times, P is 
the orbital period and Nxi is the transit number. Given all the times and transit num- 
bers, we fit for Tzp and P, and found them to be 2455307.892663±0.000082 BJDtcb, and 
1.58040487±0. 00000067 days respectively. We set the zero-point to be the first transit ob- 
served by APOSTLE (TO.Tl). In Figure E] we show the observed minus calculated (0-C) 
transit times for these transits. The plot shows that there are no significant variations from 
the expected times of transit. 

If this system did have additiona l plan ets they might induce variations in transit times 



(lAgol et al.l I2OO5I : iHolman fc Murray! l2005l ). For the special case where these additional 



planets are in resonance with GJ 1214b, the resulting TTVs could be on the order of minutes 
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and well above the precision limits of transit follow-up surveys. However such a signal does 
not seem to exist (Figure [5]). A null result however can be used to place in teresting limits 
on the mass and orbital co nfigurations of a possible undetected companion (ISteffen fc Agol 



20051 : lAgol fc Stefferu 120071 ). We will defer such an analysis until more data are available, as 



GJ 1214 remains part of APOSTLE's observing program. 



5. Absolute Stellar and Planetary Properties 



Transit and RV data when taken together allow us to constrain various planetary param- 
eters. However, many of these directly depend on the estimate of absolute stellar parameters, 
the most important being absolute stellar mass and radius. Transit lightcurves allow us to 
constrain the average stellar density (p*). We can rewrite Kepler's 3rd law to get an ex- 
pression of stellar density p-^ = 37i{a/ R^ ,)^ /(GP"^), where G is the universal gravitational 



constant (ISeager fc Mallen-Ornelas 



20031 ). The parameters, a/i?^ (the semi- major axis of 



the planet in s tellar radius units) , and P (the orbital period) can be deduced from transit 
data (see § l3.3UCarter et al.ll2008l ). To translate this measurement into an estimate of abso- 



lute mass and radius an additional constraint is required. Usually one seeks an estimate of 
the stellar mass. Once the stellar mass is obtained, getting the stellar radius is trivial, since 
= 3MJ{47tRI). 

Stellar mass is usually obtained using empirical or theoretical mass-luminosity relations 
of stars. However, given the stellar density measurement, one may also constrain mass and 
radius from the locus of points where the measured density (from transits) intersects with 
well-known stellar mass-radius relations. Mass-radius relations provide an understanding of 
the internal structures of stars. For stars in the mass range of G J 1214, empirical mass-radius 
relations are difficult to interpret due to biases in the survey sample and large uncertainties 
in the measurement of stellar parameters. Theoretical considerations of internal structure 
show th at the overall size of a star might be strongly affected by convection and magnetic 
activity (jChabrier et al.l 120071 ). Mass- luminosity relations on the other hand stem from our 
understanding of the energy production in stars. Energy production rates from nuclear fusion 
in the core are better understood and hence the luminosity of low-mass stars relate better 
with mass (jChabrier fc Baraffelll997l : lHillenbrand fc Whitdl2004j ) . In light of this, we describe 
the derivation of absolute stellar properties of GJ 1214 using mass- luminosity relations and 
discuss how the derived radius of GJ 1214 compares to well-known mass-radius relations of 
low-mass stars. Once the absolute stellar parameters are obtained, various measurements 
from RV and transit observations can be used to derive absolute parameters of GJ 1214b. 
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5.1. Stellar Mass and Radius 



Photometry is ava ilable in 8 wavebands (UBVRIJHK) for GJ 1214 ( iDawson fc Forbes 



19921 : ICutri et al.ll2003f ). We also obtained unpublished IRACl and IRAC2 flux estimates 
from Jean- Michel Desert (priv. comm.). Together these data cover the peak of GJ 1214's 
spectral energy distribution (SED). This allowed us to flt the observed SEP with spectropho- 
tometry (derived usi ng the technique descri bed in llvezic et al.l 120071 : iMaiz ApellanizI 120061 ) 
from model spectra (IHauschildt et al.lll999l ). The errors on the optical photometry were 
increased to 15% and the infrared data to 5% to compensate for inaccuracies in spectropho- 
tometry. We flt for the effective temperatures (Tg//) and logg over the model grid. In 
addition we flt for a constant, R^^/d, which is the ratio of the stellar radius to the distance 
of the star from Earth. This factor is related to the solid angle of the star and scales the 
synthetic spectra to the observed fluxes. Our best flt produced a of 10-5 given 7 degrees- 
of- freedoms (see Figure [6]). The uncertainties on the flt parameters were derived using the 
MCMC technique described in § 13. 2[ Results from the MCMC run are shown in Figure [71 
We integrated the resulting best-flt spectra (extracted at the best-flt T^jf and logg on the 
grid) over all wavelengths and scaled it with the solid-angle (related to Ri,ld) to der ive the 
observed flux {Fobs)- We then used the parallax quoted by Ivan Altena et al.l (Il995h to get 
an estimate of GJ 1214's bolometric luminosity, = 0.0028 ± O.OOO4L0. Parameter values 
and uncertainties for all flt and derived parameters are quoted in Table |H 

With an estimate of L^, theoretical or empirical m ass-luminosity relati ons can be used 
to estimate the stellar mass. We used data presented by lBaraffe et al.l ( 119981 ) for solar metal- 
licity stars at an age of 5 Gyr and found GJ 1214 to have a mass = 0.153 ± O.OIOMq. 
GJ 1214's age is assumed to lie somewh ere between 3-10 Gyr since its measured kinemat- 
ics place it in the old disk population ( iReid et al.l Il995l ). The spread in mass- luminosity 
over this age range is insigniflcant (<^ 1%) so mass and luminosity data from the 5 Gyr 
isochrone served very well. GJ 1214's metallicity has not been directly measured, but in- 
direct means indicate that [F e/H] > 0.0. For exarnple, photometric proxies indic ate that 
its [Fe/ H] = +0.03 or +0.28 Jjohnson fc AppsI [20091 ISchlaufman fc Laughlinlboioi respec- 



tively) . iRoias-Ava 



index (I Covey et al 



a et al. J2010[ ) use near-IR equivalent widths of NAI, Cal and the H2O 
2OI0I ) for M-Dwarfs with and without planets to derive a spectroscopic 
metallicity in d icator for M-Dwarfs. Using this method they report an [Fe/H] = +0.39. The 



Baraffe et al.l (Il998l ) models unfortunately do not cover this high metallicity range, so as a 
ch eck we also e s timat e the stellar mass using an empirical mass-luminosity function presented 
bvlScalo et al.l J2007h. This flt is based on masses derived from observations of binary stars 
( iHillenbrand fc Whitd 120041 ). and the sample should represent a wide range of metallicities. 
Using this flt we derived M^, = O.148M0, which is O.OO5M0 lower than the value derived 



from the theoretical relation but well within the uncertainties. ICharbonneau et al.l (120091 ) 
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use a similar technique but use a mass- l uminosity function c alibrated for mid-infrared K- 
band luminosities to determine the mass (jPelfosse et al.ll2000[ ). Their estimate of the mass 
was 0.157 ± 0.019 Mq and is consistent with our measurement. Estimating the stellar radius 
once the mass is estimated is trivial and Figure |H] demonstrates how this is done. The two 
dashed- dotted curves on the plot of stellar mass vs. radius (shown in red and blue in the color 
version), are the la contours of constant s tellar density obtained from transit observations 
in this work and ICharbonneau et al.l (120091 ) respectively. The point where the mass estimate 
intersects the measured density gives us the stellar radius, and the points of intersection 
with the la density contours determine the errors on absolute stellar mass and radius (see 
Table SD. 

Figure m shows that our estimates for the mass an d radius of GJ 1214 do not lie close to 
the e mpirical mass-radius relations for low-mass stars (jPemory et al.ll2009l : iBayless fc Orosz 
20061 ): the empirical relations are represented by the dashed lines in the figure. There is also 
disagreement between the derived values and the theoretical relations when there is no spot 
coverage. The dark shaded region on the figure shows the s pread in mass and radius due to 
a range of metallicities and ages from theoretical models (IBaraffe et al.l Il998l) and no star 



spots. Age and metallicity alone cannot account for the discrepent data. IChabrier et al. 



( 120071 ) discuss how the internal structure of low-mass stars are affected by magnetic activity 
and convection. For example, the overall sizes of stars might be affected by cool spots on the 
surface where magnetic field lines penetrate deep into the convective layer. This situation 
can reduce the efficiency of convective energy transport, ca using the star to sett le to a larger 
radius for a given lumi nosity and fiduc i al tem perature ( [Morales et al.l l2010l ). Using the 
formalism introduced by IChabrier et al.l ( l2007l ) we computed the radii for the extreme case 
when the surface of a given star is assumed to be completely covered by spots which are 
cooler than the fiducial surface temperature by 500K. Figure [8] shows the region bounded 
by this extreme as the light shaded region. The variation in radius for a given stellar mass 
is a very strong function of spot- coverage. The derived absolute parameters fall within this 
region and suggest that GJ 1214 may have cool regions on its surface. GJ 1214 would 
have to have 92% of its surface covered by spots 300K cooler than its surface to explain 
the estimated radius. For spots cooler than the surface by 500K, the large radius may be 
explained by spots covering 61% of the stellar surface area. We discuss some evidence for 
spots and magnetic activity in a subsequent section (§|6]). Addressing convection is beyond 
the scope of this paper. 

An alternate explanation for the discrepancy could be the estimate of GJ 1214's lu- 
minosity. If empirical mass-radius relations are to be believed GJ 1214's mass and radius 
would be closer to 0.23 and 0.24 Rq, based on the intersection with the density contours. 
Working backwards, this translates to a bolometric luminosity of roughly 0.0065 Lq. This 
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is a very large difference in luminosity from the current estimate. The greater luminosity 
can only be reconciled with t he observed fluxes if it w as further away by 7pc (i.e. a parallax 
that was smaller by 0.027"). Ivan Altena et al.l (Il995[ ) measure the parallax for GJ 1214 to 
be 0.0772±0.0054". The discrepent measurement would have to be a 5a systematic error 
and hence highly unlikely. 



5.2. Absolute Planetary Parameters 

With the semi- amplitude [K) from RV measurements, and the Period (P) and inclina- 
tion from transits observations, one can determine the planetary mass: K oc Mp sin i 
The planetary radius can be determined from the planet-to-star radius ratio measured from 
the transit lightcurve. From the planetary mass and radius, we derive the planet's density 
(pp), escape velocity {Vescp) and surface gravity (gp). Using the semi- major axis (a/Ri,), 
and stellar effective temperature {Tgff) we estimated the equilibrium temperature of the 
planet assuming a Bond albedo of 0.0 and 0.75. All errors were propagated assuming Gaus- 
sian uncertainties. We list our estimates in Table HI The planetary mass, radius and density 
we derive are Mp = 6.37 ± 0.87Me, Rp = 2.74tomRe, and Pp = 1.68 ± 0.23g/cm^, respec- 
tively. These data confirm GJ 1214b's status as a Super-Earth and the density measurement 
attests the presence of a massive gas envelope. 



6. Stellar Activity 



Many main-sequence stars are believed to be magnetically active and the frequency of 
active stars is known to increase with decreasing mass. Most observed variability on suc h 
stars is likely due to star spots and stellar flares (IBasri et al.ll2010t IWalkowicz et al.ll2010l ). 
The evidence for spots on the stellar surface is inferred from rotationally modulated long- 
term periodic trends in stellar lightcurves. Flares on the other hand are short-term events 
which are believed to be caused by the sudden release of energy from the reconnection of 
magnetic field lines near an active surface region. Some of the most active M-Dwarfs are 
believed to l ie at masses below the transition between partially and fully convective interiors 
(< O.35M0. iReiners fc Basrill2009l ). The increased activity is often thought to be a result of 
asymmetric magnetic field topologies for fully convective low mass stars. Although GJ 1214 
(0.153Mr7^) lies well within this mass range, it has been classified as an inactive M-Dwarf by 
Hawley et al.l (Il996l ) based on a Ha activity index. Active stars are bound to have numerous 
cool or hot regions on their surface; the existence of such regions can be inferred from the 
detection of spot or flare features in lightcurves. As we have already discussed in § O the 
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absolute size of GJ 1214 is not in agreement with our understanding of the radii of spot-free 
stars. 

We beheve that our r-band observations show evidence for a low-energy stellar flare 
on GJ 1214 and the possible detection of a cool spot on the surface. The sharp rising and 
falling trend seen in the out-of-eclipse lightcurve on UTD 2010-04-21 (Figured]) is similar 
to the fast-rise exponentia l decay (FRED) shape commonly associated with stellar flares 
(IHawley fc PettersenI Il99ll ). Panel (a) in Figure [H] shows this event in greater detail. We 
built a lightcurve model with two components: 1) a linear rise phase and 2) an exponential 
decay phase. We fit for the start time, peak time, peak flux and the e-folding time of the 
exponential phase. The best-fit lightcurve is shown as the solid gray line in Figure [9] (a). 
The Ax^ for this flare model compared to a straight line fit to the data is 112.6. The r- 
band flux of the star rose to a peak 0.8% above the quiescent level and decayed over ~3 
minutes (e-foldi ng time). Since we la cked flux-calibrated photometry, we used synthetic 
stellar spectra (IHauschildt et al.l 119991 ) to estimate the energy output by this event. We 
determined the r-band flux by integrating the synthetic sp ectra of a star with T^ff = 2949.ft" 
and l og q = 4.94 over the spectral response of the r-band (jivezic et al.ll2007l : iMaiz Apellaniz 



20061 ). Using a stellar radius of 0.21i?Q, we computed the quiescent luminosity in the r- 
band to be ~ 1.6x10^^ ergs/sec. Panel (b) shows our flare model in lu minosity units above 
the quiescent level f or GJ 1214. Following the method described in iHawley et al.l (|2003[ ) 
and iKowalski et al.l ( 120101 ) we integrated under the flare lightcurve and estimated the total 
energy output by the flare in the r-band to be ~1.8xl0^'^ ergs, see Figure IHI^c). The time it 
would take for the non-flaring star to emit this amount of energy (refered to as the equivalent 
time in the M dwarf flare community) in is 0.113 seconds. Compared to typical M-Dwarf 
flares, this event is short-lived and of much lower energy. In fact such events are likely to 
be drowned out by noise for most flare monitoring cam paigns as milli-mag p recision is not 
commonly desired when looking at the most active stars. IHawley et al.l (120031 ) reported flare 
energies between 8-58 x 10^° ergs from Johnson R-fllter observations of the active star AD 
Leonis. The activity observed on GJ 1214 is 4 orders of magnitude lower in energy than some 
of the energetic flares observed on AD Leo (also an M4.5V star). AD Le o has been identified 
as a member of the young galactic disk population (IMontes et al.l 1200 ll ). IWest et al.l (120081 ) 
have established that stellar activity decreases with age, and hence the differences in the 
activity levels of GJ 1214 and AD Leo might be purely due the differences in their ages. 

During the transit on UTD 2010-06-06, we observed a slight brightening in the lightcurve 
signal at the onset of egress (see Figure [T]). Figure [TD] shows this event in greater detail. 
The brightening could be attributed to either another flare event or the planet's temporary 
obscuration of a cool region on the stellar surface (spot-crossing). The flare model described 
above provided poor flts to this signal; the shape seen is far more symmetric than the FRED 



- 19 - 



shape of a flare. The symmetry and the fact that it occurred during transit makes it very 
hkely that we observed a spot -crossing event. We na o deled the spot-crossing signal based 
on the analytic expressions in iMandel &: Agoll (|2002[ ) for the area of intersection between 
two circles, assuming the spot was of roughly circular shape. We did not account for the 
deformation of the spot due to the curvature of the star. We fit for a parameter set very 
similar to that used for the transit model. Our variables include the spot-ingress duration, 
the spot-crossing duration, the central crossing time, the square of the radius ratio of the spot 
to the planet {R^p/Rp) and the spot temperature (T^p). To account for the fact that the cool 
region may not be a purely dark spot, we used a contrast ratio dependent on the the stellar 
and spot temperatures. So the height of the spot signal is r oughly equa l to R^^p/ R^{1— I sp / h) , 
where the contrast factor hp/h = B{Tsp, Xobs)/ B(T^,, Xobs) (!Silvall2003l ). The function B{T, A) 
is the Planck function; for this calculation the spot and star were assumed to have blackbody 
SEDs. The best fit model is shown in gray in Figure [101 The model fits the spot feature 
better than a flat-line, with the relative goodness of fit, Ax^ = 33.3. We find that a circular 
spot with a temperature ~ 2769i^ and spot-to-planet radius ratio of 0.06 fits the feature 
very well. Using the planetary radius estimated from chain 003a in Table [3] we estimate the 
radius of the spot to be ~ 0.16/?®. The duration of spot-crossing (tr.sp) can also be used to 
estimate the longitudinal extent of the cool region, Igp oc tT,spCt/PeYiod ~ 0.99i?©. The hugely 
different estimates of spot size show just how degenerate the various parameters in our model 
are. There is not enough information to reconcile degeneracies in the planet-to-spot impact 
parameter, the spot radius and the contrast ratio between the spot and star. 



7. Conclusions 



I. A Transit Model Suited for Bayesian Analysis: We show that fitting for the transit 
duration (tx) and the ingress/egress duration {tc) results in a parameter set with few mutual 
degeneracies (see Figure [2]). This condition is suited very well for MCMC methods, which 
are regularly used to determine uncertainties on parameters derived from transit lightcurves. 
Our joint analysis of multi-wavelength data using this parameter set was able to reproduce 
previous estimates of system parameters for GJ 1214b (see Table |3l chain003a). We also 
find that milli-magnitude photometry may not be sufficient to constrain limb-darkening 
parameters using transit lightcurves. We show that MCMC runs where we fit for these 
parameters were slow to converge (see Tabled]), and posterior probability distributions for 
various parameters were plagued with degeneracies (see Figure E]). Estimates of system 
parameters from these runs were generally unreliable when compared to runs where the 
limb-darkening parameters were kept fixed (see chains 'b' vs. 'a' in Tables |2] and [3]) . 
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II. Transit Timing Variations: Data gathered so far do not indicate significant varia- 
tions in the times of transit for GJ 1214b (see Figure [5]). APOSTLE will continue making 
observations of GJ 1214b and a more detailed analysis of timing data will follow in a future 
paper. 

III. System Parameters for GJ 1214: From fitting SEDs to photometry, we constrained 
GJ 1214's observed fiux and luminosity (§ |5]). The luminosity allowed us to constrain GJ 
1214's mass and since we obtained stellar density from transit lightcurves it allowed us 
to estimate GJ 1214's radius. We find the derived values of mass and radius to be in 
agreement with previous estimates, however we find GJ 1214 deviates from well-known mass- 
radius relati ons for low-mass star s (see Figure [8]). Simple calculations using the formalism 
presented in lChabrier et al.l ( 120071 ) show that GJ 1214's position on the mass-radius plot can 
be explained by the presence of cool regions on its surface. 



From RV, transit data (ICharbonneau et al.l |2009[ ) and absolute stellar properties we 
determined various properties of GJ 1214b (see Table H]). The planetary mass and radius 
(6.37±0.87 Me, 2.74;°;°^ R®) places GJ 1214b between the terrestrial and ice-giant regime 
of planets (2Me < Mp < lOM^). Its classification as a "Super-Earth" remains and the 
planetary density confi rms it is not like the rocky bodies of our solar system (see Table H]). 
Rogers fc Seagerl ( 120101 ) propose 3 scenarios for the origin of its gase ous envelope: i) primor- 
dial H /He, ii) sublimated ices (H20,C02) or iii) volcanic outgassing. iMiller-Ricci &: Fortney 
( 120101 ) propose that space-based observations of the transmission spectra of GJ 1214b's at- 
mosphere should be able to tell us how Hydrogen-rich its atmosphere is. The largest source 
of uncertainty in our estimate of planetary mass was the velocity semi-amplitude (K). Er- 
rors in the planetary radius follow from our uncertainty in measuring the absolute size of 
the star, which ultimately hinges on our luminosity estimate (see Improved precision 
on radial velocity, fiux and distance would tighten our constraints on the absolute mass and 
radius of GJ 1214b. 

IV. Evidence for Stellar Activity: The detection of a low energy stellar fiare and the 
possible transit of the planet over a star-spot (see Figures [9] and [10]) indicate that GJ 1214 is 
active. However, considering its age and comparing the fiare ener gy to fiares on the younge r 
AD Leo confirms that GJ 1214 is a quiet star for its spectral type (iHawley et al.lll996U2003l ). 
We find a fast-rise exponential decay profile fits the fiare signal (UTD 2010-04-21) quite well. 

A symmetric rise in the normalized fiux ratio during the transit on UTD 2010-06-06 
could indicate the planet occulted a star-spot on the surface of GJ 1214. The signal is 
weak, and our attempt at fitting a simplified spot model shows that spot properties are 
difficult to constrain from a single spot-crossing observation. Our results on spot-properties 
are inconclusive due to degeneracies between spot-size, planet-to-spot impact parameter and 
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spot-to-star contrast ratio in our model. Detections of this signal from successive transits 
would have confirmed it as a star-sp ot and provided intere sting constraints on the properties 
of an active stellar surface region (jPittmann et al.ll2009l ). The stellar rotation rate might 
have also been estimated with such data. 
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GJ 1214b Transits from APOSTLE 
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Fig. 1.— (a) Lightcurves of 3 transits of GJ 1214b observed by APOSTLE. The vertical 
axis shows the normalized flux ratio and the horizontal axis shows time from mid-transit 
time in days. The transit time (TO) estimated for each transit was subtracted from the time 
stamp for each point. The solid circles represent the points used in the analysis, while the 
open squares represent the data that were excluded. The excluded data were either part 
of the stellar flare (UTD 2010-04-21) or the possible spot-crossing event (UTD 2010-06-06). 
The gray line shows the transit model. Panels (b), (c) and (d) show the residuals from the 
lightcurves normalized by the transit model. The typical scatter in the normalized flux ratio 
(excluding the flare, spot and transit signal) was ~ 0.0011. 
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chainOOla Model Parameters 
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Fig. 2. — Joint probability distributions for all fitted transit 'Model' parameters fit in MCMC 
chainOOla. The model parameters were fit to APOSTLE lightcurves using the parameter 
set Oi- The numbers and units correspond those listed in Table |2l The solid-line crosshairs 
mark the location of the best-fit values (also from Table [2]). The contours mark the 1, 2, 3, 
4 and 5 sigma regions, each enclosing 68.27%, 95.45%, 99.73%, 99.994% and 99.99994% of 
the points in the distributions respectively. 
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chainOOla Derived System Parameters 
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Fig. 3. — Joint probability distributions for all 'Derived' parameters fit in MCMC chainOOla. 
The numbers and units correspond those listed in Table [2j The solid-line crosshairs mark 
the location of the best-fit values (also from Table [2]). The contours mark the 1, 2, 3, 4 and 5 
Sigma regions, each enclosing 68.27%, 95.45%, 99.73%, 99.994% and 99.99994% of the points 
in the distributions respectively. Several parameters have strong mutual correlations. 
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chainOOlb Model Parameters 
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Fig. 4. — Joint probability distributions for all transit model parameters fit in MCMC 
chainOOlb. The model parameters were fit to APOSTLE lightcurves using the parameter 
set d\. The numbers and units correspond those listed in Table El The solid- line crosshairs 
mark the location of the best-fit values (also from Table [2]). The contours mark the 1, 2, 
3, 4 and 5 sigma regions, each enclosing 68.27%, 95.45%, 99.73%, 99.994% and 99.99994% 
of the points in the distributions respectively. This set is different from chainOOla, as the 
limb-darkening parameters vIapostle and v2apostle were added to the analysis. There 
are strong correlations seen between various model parameters (e.g. D vs tj- and vs 
VIapostle)- We can also see how the non-linear minimizer converges to a low value for to 
and a high value for tr- This seems to be a result of a degeneracy with the open vIapostle 
parameter. 
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Fig. 5. — The observed minus calculated (0-C) transi t times versus transi t nurn ber for GJ 
1214b. AH transit times were converted to TCB. The ICharbonneau et al.l (12009) tira es are 
from our analysis of the joint APOSTLE and MEarth dataset. The ISada et al.l (|2010( ) data 
were converted to BJDtcb- 
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SED GJ 1214 
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Fig. 6. — The spectral energy distribution (SED) of GJ 1214 and the resulting best-fit 
spectrophotometry as described in § [51 The photometric errors on the optical (UBVRI) data 
were adjusted to 15%, while those for the infrared data (2MASS and Warm Spitzer) were 
adjusted to 5%. 
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SED Fitting Parameters 
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Fig. 7. — Joint probability distributions for the three parameters used for SED fitting. The 
numbers and units correspond to those hsted in Table |H The solid-line crosshairs mark the 
location of the best-fit values (also from Table H]). The contours mark the 1, 2, 3, 4 and 5 
Sigma regions, each enclosing 68.27%, 95.45%, 99.73%, 99.994% and 99.99994% of the points 
in the distributions respectively. 
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Stellar Mass vs. Radius 
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Fig. 8. — Stellar mass (M^) vs. stellar radius (R*), in solar units. The linestyles are 
matched to their corresponding references in the legend on the top-left corner of the figure 
and the legend on the bottom right matches the data points on the plot. The two dashed 



Demorv et al. 


2009; 


Bavless & Orosz 


2006) 



and ra dius. The estimate of mass and radius presented in this work and ICharbonneau et al. 



(120091 ) are also marked on the plot. The two dashed-dotted curves, shown in red and blue 
in the color version, are the l a contours of const a nt ste llar density obtained from transit 
observations in this work and ICharbonneau et al.l (120091 ) respectively. The shaded regions 
represent the spread in stellar mass and radius taking into account various theoretical con- 
siderations. The darker region represents the spread over age (3-10 Gyr) and metallicity 
(-0.5-0.0) for stars without spots. T he lighter regio n shows this spread when spot-coverage 
is introduced using the formalism of IChabrier et al.l (120071 ). The outer limit (leftmost edge) 
of the light region represents the extreme case of 100% area coverage of spots which are 500K 



coole r than the fiducial the surface temperature of the star for a given mass (IMorales et al. 
2010h . 
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Fig. 9. — (a) Shows the flare event observed on UTD 2010-04-21. The data are the same 
as those shown in Figured] (b), only magnified. The gray line is the best-fit FRED model. 
Even though only the points after mid-transit are shown, the fit was made using all lightcurve 
points (in Figured] (b)). Panel (b) in the above figure shows the r-band luminosity as a 
function of time and panel (c) shows the total energy output by the fiare above the quiescent 
r-band level as a function of time. 
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Possible Spot Crossing During GJ1214b Transit 
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Fig. 10. — Possible spot-crossing event from the UTD 2010-06-06 transit of GJ 1214b. The 
transit signal was removed by normalizing the lightcurve with the best-fit transit model. 
The vertical dashed lines approximately mark start and end of the transit. The gray line 
shows a fit using a simplified spot model. 



Table 1: MCMC Analysis 



Chain Name 


Data Set 


Model Pars 




Chain Length 


Corr Length 


Eff Length 




DOF 


001a 


APOSTLE 


01 


6 


998,301 


11 


90,755 


802 


808 


001b 


APOSTLE 


01 


8 


1,193,465 


354 


3,371 


796 


806 


002a 


MEarth + FLW01.2m 


92 


9 


998,308 


66 


15,126 


1279 


1294 


002b 


MEarth + FLW01.2m 


92 


13 


1,193,569 


939 


1,271 


1249 


1290 


003a 


APOSTLE + MEarth + FLW01.2m 


92 


12 


997,202 


84 


11,871 


2067 


2105 


003b 


APOSTLE + MEarth + FLW01.2m 


92 


18 


2,398,779 


3,299 


727 


1947 


2099 


004a 


MEarth 


9i 


7 


996,616 


16 


62,289 


804 


813 


004b 


MEarth 


9i 


9 


1,195,033 


240 


4,979 


796 


811 


005a 


FLW01.2m 


9i 


5 


997,093 


8 


124,637 


475 


478 


005b 


FLW01.2m 


9i 


7 


1,198,883 


148 


8,101 


471 


476 



Table 2: Parameter set di 



Parameter/Chain 


001a 


001b 


004a 


004b 


005a 


OOSb 


Charbonneau et al. (2009)* 


Units 


Model 


D 


0.0135 ± 0.0002 


0130+0-OUlu 


0.0141 ± 0.0001 




0.0145 ± 0.0003 


n ni OC + 0.UUU3 

0-0136-0.0005 


0.0135 ± 0.0002 




tT 
tG 


0.0326 ± 0.0002 


0329+0-0001 
n nnis+O-0014 


0.0311 ± 0.0001 
0037+0-0009 


0315 + 0-0003 
0040+0-0006 

U.UU4U_Q QQQ^ 


0.0314 ± 0.0002 

0051+0-0009 

U.UUai_(j ggg2 


0319+0-0006 

U.Uai»_g ggg^ 

0045+0-0010 
0-00*^-0.0004 


0.0321 ± 0.1529 
0.0043 ± 0.0203 


day 
day 


^'^APOSTLE 


(0.908) 


1 00+0-06 














V^APOSTLE 


(0.305) 


42+0.59 
"■1^-0.07 




















(0.145) 


0.28 ± 0.12 






(0.145) 




"^'■■'Ea^th 






(0.639) 


1 24+0-^^2 






(0.639) 




vIfL W01.2T7^ 










(0.404) 


O-°l-0.26 


(0.404) 




v2fL W01.2m 










(-0.289) 


30+0-03 
0--50_Q 57 


(-0.289) 




TO.Tl 2455307+ 


0.892740 ± 0.000275 


0.892769 ± 0.000276 












B.JDrcB 


T0.T2 2455353+ 


0.724696 ± 0.000317 


0.724738 ± 0.000317 












B.JDrcB 


T0.T3 2455383+ 


0.752366 ± 0.000265 


0.752306 ± 0.000264 










00 


B.JDrcB 


T0.T4 2454964+ 






0.944783 ± 0.000794 


0.944955 ± 0.000792 






0.945142 ± 0.000403 ^ 


BJDrcB 


T0.T5 2454980+ 






0.748909 ± 0.000272 


0.748904 ± 0.000271 






0.748903 ± 0.000090 


BJDtcb 


T0.T6 2454983+ 






0.909659 ± 0.000233 


0.909675 ± 0.000233 






0.909693 ± 0.000090 


bJDtcb 


T0.T7 2454999+ 






0.713612 ± 0.000260 


0.713654 ± 0.000260 






0.713634 ± 0.000126 


BJDtcb 


T0.T8 2454980+ 










0.748974 ± 0.000467 


0.748886 ± 0.000469 




BJDtcb 


T0.T9 2454983+ 










0.909754 ± 0.000444 


0.909699 ± 0.000446 




BJDrcB 


Derived 


Rp/ R* 


0.i2J4_(, pp^-. 


1143+0-0134 
U.J-140_Q QQ33 


„ „_+U.UU42 
O-llS'^0.0006 


-,-is2 + 0-0U48 

0.11»2_g ggj^j; 


1221+0-00^0 

U.lZZl_g ggj^g 


1205+0-006S 

U.i^UO_g ggg^ 


0.1162 ± 0.0007 




P(1.5804 days+) 


0.79 ± 0.25 


0.65 ± 0.25 


0.10 ± 1.38 


-0.40 ± 1.39 


-0.86 ± 9.97 


0.52 ± 9.91 


-0.65 ± 1.01 


sec 


b 

a/fl« 
i 

V 

P* 


41 + 0-03 
14 06+''-''^ 

55.901=?^ 


11 + 0-35 

15.18l5'.26 
89.5915:19 
60 34+1-07 
26 49+1-3^ 


04+0-3S' 

16.191°;- 

89.87l?;f3 
64.3510-85 
32 13+"^ 26 


25+0-1* 

i5.47lg;- 
89.08l°o;l^ 

61.52lO;7^, 
28.06li;0/g 


o.5olS;i°3 
I3.90i;;f3 

87 94+0-17 
°'-^*-0.61 

55 25+1-02 
20-3414.19 


„ ,„+0.14 

0.3»_o 10 
1/1 R1+0-43 

14.61_o gg 

88.5210;** 
58.0711;^! 
23.6il2;i= 


35+0-00 

14.66 ± 0.41 
88.62l0;3| 

58.28 ± 1.63 
23.90 ± 2.10 


deg 
day-1 
g/cm3 



* For those parameters not explicitly listed by ICharbonneau et al.| j2009l l (e.g. t'j' and tQ) wc computed them using the expressions in ICarter et aT] j2008l ^ and propagated the errors 
assuming they had Gaussian distributions. 



Table 3: Parameter set 62 



Parameter/Chain 002a 002b 003a 003b Charbonneau et al. (2009) ' Units 



Model 



tT 
tG 

'^'^ APOSTLE 


0.0145 ± 0.0002 
0310+''-'""'! 
0045+'' ™''= 

U.UU40_Q QQQ3 


0137+"-'""" 

0316+° °°°^ 
"-"•^-'°-0.0004 


0143+"-'""*^ 
0.0316 ± 0.0001 
0045+0-'""'^ 
(0.908) 


„ „,,,+O.UU(lti 
"-"^■^■^-0.0001 
0321+0 0002 

0.0037 ± 0.0002 
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'-^-0.08 


0.0135 ± 0.0002 
0.0321 ± 0.1529 
0.0043 ± 0.0203 


day 
day 


APOSTLE 






(0.305) 


1 01+0-20 








(0.145) 


0.38l°-0^ 


(0.145) 


53+0-11 


(0.145) 




vli^JLW01.2m 


(0.639) 
(0.404) 




(0.639) 
(0.404) 


47+0-6S 
-0.28 

61+0-13 
"-"4-0.20 


(0.639) 
(0.404) 




^^F LW01.2'm. 


(-0.289) 


0.55 ± 0.58 


(-0.289) 


0.24 ± 0.51 


(-0.289) 




TO.Tl 2455307+ 






0.892642 ± 0.000271 


0.892689 ± 0.000263 






T0.T2 2455353+ 






0.724728 ± 0.000307 


0.724652 ± 0.000311 




BJDtcb 


T0.T3 2455383+ 






0.752332 ± 0.000260 


0.752334 ± 0.000264 




BJDTCiJ 


T0.T4 2454964+ 


0.944943 ± 0.000794 


0.944940 ± 0.000789 


0.944962 ± 0.000799 


0.944935 ± 0.000788 


0.945142 ± 0.000403 


BJDtcb 


T0.T5 2454980+ 


0.748921 ± 0.000266 


0.748938 ± 0.000267 


0.748905 ± 0.000263 


0.748976 ± 0.000264 


0.748903 ± 0.000090 


BJDtcb 


T0.T6 2454983+ 


0.909665 ± 0.000229 


0.909650 ± 0.000230 


0.909674 ± 0.000228 


0.909689 ± 0.000228 


0.909693 ± 0.000090 


BJDtcs 


T0.T7 2454999+ 


0.713648 ± 0.000256 


0.713675 ± 0.000257 


0.713625 ± 0.000253 


0.713690 ± 0.000253 


0.713634 ± 0.000126 


BJDrciS 


T0.T8 2454980+ 


0.748938 ± 0.000380 


0.748897 ± 0.000414 


0.748924 ± 0.000390 


0.748942 ± 0.000417 






T0.T9 2454983+ 


0.909734 ± 0.000373 


0.909715 ± 0.000400 


0.909736 ± 0.000381 


0.909686 ± 0.000401 




BJDrcfs 


Derived 


P(1.5804 days + ) 


U.J-ZUO_Q QQQg 

-0.39 ± 1.38 


U.1J-/J-_Q QQQg 

-0.27 ± 1.42 


1195+"-'"""* 
0.47 ± 0.04 


1152+0-0020 
"--4-4^^-0.0003 

0.46 ± 0.04 


0.1162 ± 0.0007 
-0.65 ± 1.01 


sec 


b 

a/ii* 

V 

P* 


14.70i«:« 
24-09i^:™ 


0.3210-.J3 

15.08lS:f, 
88.79lg;l| 
59 94+1 10 
25-97tlJ| 


14 54+0-10 
88.39+0-1* 
57.81+0-^^ 
23.29t;o/5 


09+0-2* 
"■"^-0.04 

1 e cn+0.08 
I5.6OI0 ,^2 

89 67+0-15 
62.02+0 1^ 
28.76lo;« 


35+0-06 
"-•^"-0.08 

14.66 ± 0.41 
88.62+2-JI 
58.28 ± 1.63 
23.90 ± 2.10 


dcg 
day-1 
g/cm^ 



* For those parameters not explicitly listed by [Charbonneau et al.| j2009| ) (e.g. tj^ and t^^) wc computed them using the expressions in [Carter etlHI i2008l ) and propagated the errors 
assuming they had Gaussian distributions. 



Table 4: Propert: 


ies of GJ 1214 and GJ 1214b 




Parameter 


Value 


Units 


Fit SED Parameters 


Teff 


2949t^^ 


K 


^ogg 


4 04+0-22 


log cm/s^ 




0.0155 ± 0.0003 


RqIPC 


Derived Stellar Parameters 


Fobs 


5.23 ± 0.13 


10 ^^ergs/ s/cm^ 




0.0028 ± 0.0004 






U.iOO_Q QQ9 


Me 




210+°-°°^ 

U.ZiU_Q QQ4 


Rq 


Derived Planetary Parameters 




6.37 ± 0.87 


Me 


Rp 


9 7/1+0.06 
^- '^-0.05 


R® 


Pp 


1.68 ± 0.23 


g/cm^ 


9p 


8.24 ± 1.09 






12.03 ± 0.80 


km/s 



T^q (Bond albedo = 0) 5471^ K 
Teg (Bond albedo = 0.75) 3871^ K 



